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Optical pulses at THz and mid-infrared frequencies tuned to specific vibrational resonances modulate the 
lattice along chosen normal mode coordinates. In this way, solids can be switched between competing 
electronic phases and new states are created. Here, we use vibrational modulation to make electronic 
interactions (Hubbard- U) in Mott-insulator time dependent. Mid-infrared optical pulses excite localized 
molecular vibrations in ET-F2TCNQ, a prototypical one-dimensional Mott-insulator. A broadband 
ultrafast probe interrogates the resulting optical spectrum between THz and visible frequencies. A 
red-shifted charge-transfer resonance is observed, consistent with a time-averaged reduction of the 
electronic correlation strength U. Secondly, a sideband manifold inside of the Mott-gap appears, resulting 
from a periodically modulated U. The response is compared to computations based on a 
quantum-modulated dynamic Hubbard model. Heuristic fitting suggests asymmetric holon-doublon 
coupling to the molecules and that electron double-occupancies strongly squeeze the vibrational mode. 

A Mott insulator is a solid with fractionally filled electronic bands in which charge carriers do not conduct 
because of their mutual repulsion. The resulting properties^ are believed to be captured by Hubbard 
models^'^, where only the screened Coulomb interaction, i.e. the on-site Hubbard- (7, is retained. This 
approach has been widely applied to polymers, transition metal oxides, fullerenes and cuprates^"^. Moreover, the 
Hubbard model can effectively describe optically excited Mott insulators where hot quasi-particles are injected 
across the gap. The quantum dynamics following such "photo-doping" has been extensively studied in the solid 
state^°-^l 

Non-equilibrium Hubbard physics is routinely explored in 'artificial solids' realized by ultra-cold atoms 
trapped in optical lattices^'*"^^. However, a powerful feature of these systems is the possibility to accurately and 
dynamically control the strength of interactions through the Feshbach resonances, and the bandwidth via the 
lattice depth^^. In contrast, the ability to control the parameters of 'real' materials is far more restricted. While the 
bandwidth can be statically controlled across a limited parameter space, by using pressure^^ or chemical sub- 
stitutions, the Hubbard- U cannot easily be modified since it is dependent on the atomic physics of each site. 

A promising pathway to such control are optical pulses at THz and mid-infrared frequencies, which can be 
tuned to specific vibrational resonances and modulate the lattice along a chosen normal mode coordinate. Since 
the photon energies are far below those of the correlation gap, this technique has already been used to switch 
solids between competing electronic phases^°'^\ and to create new states not found at equilibrium^^. Here, we 
investigate how vibrational modulation influences the electronic interactions in a Mott insulator. Crucially, the 
leading-order effect of exciting these localised modes is to dynamically modulate the on-site electronic wave- 
function, thereby inducing temporal oscillations and modifying the time-averaged Hubbard- (7 interaction 
strength. 

The response of a strongly correlated electronic system to such vibrational modulation is not well understood, 
yet the dynamical states created have new properties and exhibit less dissipation than in the case of quasi-particles 
excitations, opening up the potential for optical devices with superior thermal management. 

Results 

We investigate at room-temperature an organic Mott insulator, ^/s(ethylendithyo)-tetrathiafulvalene-difluoro- 
tetracyano-quinodimethane (ET-F2TCNQ), which is a prototypical system where Hubbard physics accurately 
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describes the electronic transport between individual molecular 
sites^^. The static optical properties of ET-F2TCNQ are shown in 
figure 1 for light polarized along the chain of ET molecules. A 
charge-transfer resonance, corresponding to excitations of the type 
(ET^, ET^) (ET^^, ET°) is observed at photon energies —700 meV 
(5500 cm"^), reflecting the existence of a correlation gap. 

In the time resolved experiments reported in this work, we mea- 
sure the broadband optical properties as one localized intramolecular 
mode is driven with mid-infrared femtosecond pulses. These were 
generated by difference-frequency mixing of two near-infrared 
pulses, obtained by parametric down -conversion of a single 800 
nm wavelength femtosecond pulse. The excitation pulses were polar- 
ized perpendicular to the chains, driving an infrared active molecular 
vibration close to —10 |im wavelength (Q = 1000 cm"^)^^. This 
10 jim mode was excited with field strengths up to 10 MV/cm, thus 
strongly deforming the molecular oscillators and local molecular 
orbitals (see supplementary information). 

The tuneable output of the second optical parametric amplifier 
was used to probe the reflectivity changes along the ET molecule 
chains in the mid (1800-3000 cm"^) and near infrared (4000- 
7000 cm"^). The THz response was probed between 25 and 
85 cm"^ using single-cycle pulses generated by optical rectification 
in ZnTe. 

When excited by the mid-infrared pump pulses, a shift of the 
charge transfer band to lower frequency was observed, along with 
new distinct peaks inside the gap (figure 2a). In this driven state, no 
reflectivity change was detected for frequencies above the charge 



transfer resonance, and no metallic response was detected in the 
THz range, where the reflectivity remained low and the phonon 
resonances unscreened (see figure 2b). The lack of a low frequency 
response clearly sets this modulated state apart from the photo - 
induced formation of metallic state by above gap optical excitation, 
for which a Drude-like metallic response is typically observed^^. 
Secondly, the rapid relaxation back to the ground state makes this 
experiment different from the case of the vibrationally driven phase 
transitions, where THz frequency excitation can lead to switching 
between competing orders^^'^^. As the excited mode relaxes, the 
vibrational oscillations are reduced and the effects vanish as seen 
in figure 2a. By 0.5 ps after the excitation, the peaks inside the gap 
are strongly reduced and the charge transfer band shifts back to its 
equilibrium position. Relaxation back to the ground state occurred 
with a double-exponential decay (xi = 230 fs and T2 = 4.5 ps). The 
first fast time constant we interpret as the direct relaxation of the 
system in the vibrational excited state. The second long time constant 
we assign to a general thermalization of the hot mode with the whole 
crystal lattice that also involves low frequency phonon modes at THz 
frequencies resulting in the ps-time constant. The response is also 
mode selective. When the excitation wavelength was tuned away 
from the frequency of this vibration, the mid-infrared resonances 
disappeared. 

Quantitative analysis of the transient optical properties starts from 
fitting the broadband reflectivity spectra of figure 2b with a Drude- 
Lorentz Model and transforming it into optical conductivity^^. The 
conductivity lineouts are shown in figure 4c. The red- shift of the 
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Figure 1 | (a) Crystal structure of the ET-F2TCNQ. The ET molecules form a one-dimensional, half filled Mott-insulating chain along the a axis, (b) Black 
curves: equilibrium reflectivity (left panel) and optical conductivity (right panel) of ET-F2TCNQ for light polarized along the a-axis (ET chains). The 
optical conductivity is extracted using a Drude Lorentz fit of the reflectivity. The charge transfer peak in the optical conductivity is fitted to a 10 site 
numerical solution of the Hubbard model (red) based on the unequal time current- current correlation function (see methods). The best fit is obtained 
with hopping t ~ 40 meV, onsite interaction 17—850 meV and nearest-neighbour repulsion y ~ 120 meV. These values are in agreement with reference 
12 and the size of the optical gap in reference 23. A quantitative match with the experimental conductivity is achieved by adding a Lorentzian oscillator 
centred at 8000 cm~\ which takes into account contributions by higher-lying transitions (blue dashed curve). 
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Figure 2 | (a) Time and frequency dependent reflectivity after selective modulation of the infrared active molecular mode at 10 |^m wavelength, (b) 
Frequency dependent reflectivity at equilibrium (black) and at the peak of the modulating infrared field (red). Full circles indicate experimental data, 
whilst the dashed line is a Lorentz fit to the data to extract the optical conductivity. In the low- frequency range (below 100 cm"^), the full lines indicate the 
equilibrium (black) and transient (red) reflectivity measured with single-cycle THz-pulses. The two peaks are phonon modes of the molecular crystal. 



charge transfer band occurs from its equilibrium position at co/c ~ 
5500 cm"^ toward 5000 cm~\ and a new band is observed approxi- 
mately at 4200 cm~\ Additionally, a mid-gap resonance and a 
weaker peak appeared at —3000 cm"^ and —2000 cm~\ respectively. 

To analyse these data, we start from a description of the static 
optical properties of this compound using the extended Hubbard 
Hamiltonian, 



HhuB '- 



>(4 



+ h.c. \ + Ui:ih^hii + Vi:ihihi+i (1) 



In this expression cj^ and Qa are creation and annihilation operators 
for an electron at site i with spin cr, hia is its corresponding number 
operator and hi = hii -\- . The optical response of the system at 
high photon energies is probed by the creation of Hubbard excitons^^ 
(figure 3a), composed of a doubly occupied "doublon" site with 
repulsive energy (7, bound to an neighbouring empty "holon" site 
by a — y attractive energy. The static optical spectrum features a delta 
peak located at coct = U— V— 4f-/V, corresponding to the Hubbard 
exciton, and a broad continuum centered about (7, with a bandwidth 
of 8t, corresponding to unbound holon-doublon excitations. 
Figure lb shows that the measured spectrum is well fitted by the 
extended Hubbard model optical conductivity once a background 
contribution from a high frequency oscillator (dashed blue curve) is 
added. 

To describe the vibrationally- driven state, we consider the model 
shown in figure 3b. To each molecular site we explicitly introduce a 
harmonic oscillator describing the selected molecular mode that 
gives a locally vibrating Hubbard modeP° 



Hlvh = HhuB -\-Hi-\- T^iHi . 



(2) 



Here X^H^ is the sum of local harmonic oscillator Hamiltonians each 
with angular frequency Q, displacement qi and ground state size ao- 
The term Hi takes into account the coupling between the charge and 
driven mode, and is expected to depend on the electronic configura- 
tion as hif{qi)-\-h^hiig{qi), where /and g are general functions. 
Typical electron-phonon interactions within the Hoist ein^^ model 
would only retain a linear coupling to the charge density as hiqi. 
However, a Holstein model is not sufficient to describe the "sloshing" 
motion of the infrared active distortion (see supplementary informa- 
tion) since by symmetry linear couplings vanish. Instead we find (see 
Methods) that the vibration has a dipolar coupling to the charge 
sector which to leading order gives an interaction 



Hj = Ti{hiHi-diDi)ql (3) 

where (l — h^) (l — hii) and D£ = h^hii locai holon 

and doublon projectors, with respective positive couplings di and 

In addition to Hj coupling to the vibration quadratically, crucially 
it also depends on the double occupancy as h^fniifi making Hlvh a 
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Figure 3 | (a) Each charge transfer excitation creates an empty site (holon) 
neighbouring a doubly occupied site (doublon). The holon-doublon pair 
involves a repulsive energy U and a binding energy of V, totalling U — V. 
(b) On each site electrons in the valence orbital of the FT molecule are 
coupled to a harmonic oscillator explicitly modelling the selected 
vibrational mode. For the 10 |am mode the coupling to the electrons is 
expected to "stiffen" the oscillator when one electron is removed (holon), 
and to "slacken" when one electron is added (doublon). 
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Figure 4 | (a) Schematic representation of the vibrational transitions assuming for simpKcity = Q. The double-occupied site has a significantly 
renormalized vibration, which results in new transitions in the optical conductivity. The p„ are the occupations of the driven lower vibrational ladder, 
(b) Theoretical optical conductivity at equilibrium (as in figure lb) upon coherent excitation of the local vibrational mode. The green and red curves 
represent two different driving amplitudes. The various peaks are labelled according to the transitions of (a). The dashed curve represents the case in 
which the asymmetry between holons and doublons are inverted, i.e. in where the holons are stiffened by a factor of 4 and the doublons remain almost 
unpertubed. (c) Optical conductivity calculated via Drude Lorentz fit to the reflectivity spectra of figure 2b, at two selected pump-probe delays. The 
background oscillator of figure lb has been subtracted. 



form of dynamic Hubbard modeP^. This type of model, which is 
natural extension of the Holstein-Hubbard physics, was originally 
introduced by Hirsch^^ to energetically account for reductions in the 
Hubbard (7 arising from two electrons occupying the same site. As we 
shall show the similarity of this effect to the vibrational driving here 
necessitates the use of a dynamic Hubbard model in order to capture 
the experimental findings. 

Discussion 

To analyse the experimental data of figure 2 we computed the real 
part of the optical conductivity o-(co) of Hlyh via the unequal time 
current- current correlation function^^'^^ Xjji^)- The narrow hopping 
bandwidth of ET-F2TCNQ, as well as the strong binding of holon- 
doublon pairs, justifies g(w) being calculated analytically in the 
atomic limit where electron hopping is neglected (see Methods). 
The key effect of vibrational excitation are revealed by considering 
the infinitely heavy oscillator limit, where its classical displacement 
varies harmonically in time t as <^ = Qcos(Qt), with amplitude Q and 
Q = 1000 cm"^ the frequency of the driven molecular vibration. 
Thus, the onsite interaction matrix element U is modulated by (h 
— d)q^ and results in two crucial features. First, the cf oc Q^[l + 
cos(2I2t)] dependence predicts a shift of the charge transfer res- 
onance, by an amount that depends on the amplitude of the driven 
mode Q. The charge transfer resonance is expected to red-shift if the 
doublon coupling exceeds that of the holon, and therefore the time- 
averaged (7 is effectively reduced. A blue-shift occurs for the opposite 
case. Second, the classical frequency modulation generates sidebands 



at multiples of ±2 Q on each side of the shifted charge transfer 
resonance. That a shift in the charge transfer resonance is seen in 
figure 2a, along with the emergence of a mid-gap peak at 2 Q below, 
is strong evidence of the quadratic modulation of the interaction U 
consistent with the dynamic Hubbard model. 

A number of experimental features are not captured by the simple 
classical modulation picture. In particular a sideband is only seen to 
the red of the charge transfer resonance and additional sub-structure 
is present. Instead, to better reproduce these experimental line shapes 
with a dynamic Hubbard model implies a strong holon-doublon 
asymmetry along with a quantum treatment of a finite mass oscil- 
lator (see Methods). As shown in figure 3b, this strongly coupled 
limit results in back-action of the electronic configuration on the 
vibrational mode. Specifically, Hj describes a "stiffening" of the oscil- 
lator on the holon site, so Q/j > Q, and a "slackening" on the doublon 
site, so < Q^^'^^. Such abrupt frequency changes are well known to 
cause squeezing of the oscillators dynamics^^'^^. 

Figure 4b and 4c show the comparison with the optical conduc- 
tivity extracted from the reflectivity data and a calculation of g{(jj) 
with strong coupling aimed to reproduce the overall structure of the 
experimental conductivity data. We obtain a driving strength Q/ao ~ 
2, a doublon oscillator that suffers a significant frequency reduction 
to ~ 0.26 Q, and a holon's Q/j ~ 1.10 Q that is only marginally 
increased. Therefore the additional sub -structure seen in figure 4b 
can be attributed to the "slackening" of the doublon oscillator as 
schematically shown in figure 4a. Crucially, the reduced spacing 
between energy levels of the doublon oscillator causes the transition 



SCIENTIFIC REPORTS | 4 : 3823 | DOI: 1 0.1 038/srep03823 



4 



(a) (b) 




Time delay (ps) Wavenumbers (cm'^) 

Figure 5 | (a) Time and frequency dependent reflectivity after modulation at 6 [im wavelength, (b) Frequency dependent reflectivity at equilibrium 
(black) and at the peak of the modulating infrared field (red). Full circles indicate experimental data, whilst the dashed line is a Lorentz fit to 
the data. 



frequencies to move into the gap, an effect that is significant only if 
the vibrational mode is appreciably populated as found from the 
model parameters to fit the data. This finding is in agreement with 
the strong experimental excitation of the molecular oscillator. The 
remnants of the classical sidebands are now located only at low 
frequencies and are split into multiples of Q — Q^. We note that 
the 10% modulation of the holon frequency is in agreement with the 
expectations^^'^^ for a thermal charge ordered molecular state. 
However, the substantial reduction, by a factor 4, of the frequency 
for the doublon site goes far beyond any equilibrium predictions. 
Nonetheless, the greater effect is consistent with the charge distri- 
bution of double occupancies lying predominantly near the oscil- 
lating C=C bond^^'^^ 

By way of comparison, we also report the response when tuning 
the pump wavelength to 6 |im, near a weakly infrared- active sym- 
metric mode along the chains^°. This vibrational mode is expected to 
have a linear coupling. As shown in figure 5, only a reduction in 
spectral weight at the charge transfer resonance was observed, with- 
out any shift or significant response at other wavelengths. As 
expected, this indicates a dependence of the optical response on 
the coupling strength, driving and symmetry of the selected vibra- 
tional mode. 

This experiment demonstrates how for this compound selective 
vibrational excitation of intramolecular modes can, to first approxi- 
mation, cause a dynamical modulation of the on-site Hubbard- (7 
interaction. Such control is not easily accomplished by other means 
and the method is likely to be applicable to a broader class of organic 
materials. More refined analysis of the observed lineshape carries 
important information about the microscopic interactions of the 
material, such as holon-doublon asymmetry. Most importantly, 
the response to such strongly driven vibrational modes reveals the 
strength of their coupling to the electronic structure. 

More generally, the selective modulation of one degree of freedom, 
when combined with probing of the electronic spectrum, raises the 
tantalising prospect of experimentally deconstructing the Hubbard 
Hamiltonian, by exposing one specific coupling that would otherwise 
have a vanishingly small contribution to the equilibrium properties. 
The deconstruction obtained by such quantum modulation spectro- 
scopy should be applicable to more than one class of modes, includ- 
ing charge, magnetic or orbital excitations. 

Methods 

Coupling of local modes to electrons. In general a spin-independent coupling 
between a local harmonic oscillator and the electronic configuration of a site can be 
written in terms of hi and h^hii only, as the number operator hi is fermionic. We 
assume that this coupling can be expressed as hif{qi) + hi^^ hiig{qi), where/(^£) and 
g{qi) are two functions of the local mode coordinate that are not known a priori. By 
expanding the functions /and g into a Taylor series we obtain 



Hi = i:ihi{Aiqi+A2q]+ ■ ■ ■) +hi^hii{Biqi + B2q]+ ■ ■ ■) (4) 

with coupling constants A,, 5, which are constrained by the symmetry of the 
molecular modes'^\ An antisymmetric infrared vibration, such as the 10-)im mode in 
the experiment, is akin to an oscillating dipole. Within the Born-Oppenheimer 
approximation this perturbation causes admixing of the valence orbital with higher- 
lying excited states of differing parity and induces an energy shift that is an even 
function of qi, meaning e.g. that A i = 0. This already precludes that such a vibration 
can be described by a conventional Holstein-Hubbard^^ type interaction. In contrast, 
a symmetric Raman vibration, such as the mode at 6 |j,m, is instead captured by an 
oscillating quadrupole, admixing higher-lying states of the same parity and causing a 
linear energy shift with non-zero Ai. 

The second term in the expression for H;, which includes coupling to the 
double occupancy, is determined by computing the Coulomb repulsion arising 
from both electrons occupying the admixed vibrational orbital | Wq) via U{q) oc 

q\i{^ q\2^n'\^ q) q) 2 ^hcrc v 12 = \h—h \ is the distance between the 
electrons. Retaining these terms is crucial to properly describing how vibrations 
modif)^ the on-site Hubbard U electronic interactions. 

For the infrared vibration the differing parity of the states in the admixture causes 
U{q) to again vary, to lowest order, quadratically with displacement q (implying 
Bi = 0), while the Raman vibration retains a linear dependence. It is found that < 0 
for the infrared mode because the admixed vibrational orbital spatially expands for 
any non-zero displacement in this simple model. We finally re-arrange terms to 
isolate the holon = (l — n^j) (l — ^^|) and doublon = h^hii couplings and to 
give the expression for Hi in the main text. 

Optical conductivity of a vibrating Hubbard chain. Since U':^kBT':^f' J U the 

system does not display any long-range Neel spin-order and does not contain 
doublons in its thermal state. It is thus well approximated by a half-filled completely 
spin-mixed state 

Ps=l,E^\^){-^\ (5) 

where \a) is a spin configuration state of the electronic Mott insulator. Given that 
ksTKHQ higher vibrational states of the localized molecular oscillators are essentially 
unoccupied at room temperature. To model the vibrationally driven state we assume 
that each oscillator is prepared instantaneously in a coherent state |a). The phase of 
the oscillator on each site is expected to be identical across the sample, however it is 
not controlled shot-to- shot in the experiment. As such a global averaging of the phase 
of a is made. We further assume that the relative phase of the oscillators where the 
holon and doublon becomes incoherent upon optical excitation. Since in the atomic 
limit these are the only relevant oscillators in the problem, this is equivalent to 
assuming they are both in a phase- averaged coherent state p^. 

To compute the optical conductivity we start from the unequal time current- 
current correlation function which is defined as XjjM = ^M^^[JMJWp] where 
0(t) is the Heaviside function and /(t) = exp (/Hi /h)} exp ( — iHz / H) is the 
Heisenberg picture current operator. The regular finite frequency optical conduc- 
tivity then follows as 

ai(a;>0)ociRe{xy(co>0)-Zy(co>0)} (6) 

where Xjji^^) ^nd Xjj{(^) are the Fourier transforms of Xjji^:) and its complex conjugate, 
respectively. 
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In the atomic limit we obtain the current-current correlation function^^ 

X^-^(a»)oc| I GQ{cQ")Gh{cQ' — CQ")Gd{co — co')dcQ'dcQ" (7) 
as a convolution of the Hubbard shift 

Go(a;)oc r 0(T)e-^'(^-^-^™)^/^^iT (8) 

J — CO 

and the onsite vibrational doublon and holon correlation functions 

Gh4H ^ ^n,mpn \ {n\s{^h,d)\m) \'s(^w-(^m^ 0 Qm + + ^) ^) 

where S[^h4) is the squeezing operator with c^^^ = log ^/OhJ/Q and pn are the 
diagonal elements of p^. The matrix elements \{n\S(^^h,d)\^)f' weighting the 5 
functions are Franck- Condon factors and are well known analytically*^ '^^. A similar 
analysis can proceed for the Raman mode. 

Atomic limit and hopping. The optical conductivity is calculated here in the so- 
called zero -bandwidth atomic limit, where the ratio Gi{oj)lf can be obtained exactly. 
This approach is justified by the parameters of ET-F2TCNQ, which indicate that an 
adjacent holon-doublon pair will be bound. Thus, once optically excited the 
contributions of the oscillators at the initial holon and doublon locations will be 
dominant. Moreover, in the dynamic Hubbard model Hlvh> the quadratic electron- 
oscillator coupling causes a suppression of t, analogous to polaronic effects^*'. A 
further suppression of the coherent hopping processes occurs once the vibrational 
modes are driven, similar to the reduction with temperature seen for polarons. For the 
strong coupling infrared case we estimate that the hopping of holons drops to around 
0.2 t, while for doublons it is just 0.05 t. Also, because the oscillator coupling here is to 
a local molecular mode on each site, as opposed to a collective bath of lattice phonons, 
vibrational excitation disorders hopping through the chain, further inhibiting holon- 
doublon motion. Thus, the effects of a finite t are to broaden the dominant 
contribution to the optical conductivity already captured by the on-site vibrational 
dynamics. To account for these mechanisms, as well as the spectral limitations of the 
measurement itself, we introduced an artificial broadening of 0.5 tio the results 
presented. 
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